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Abstract 


In this paper, we explain a new Iterative Method - Fixed Point and develop its convergence theory 
for finding approximate solutions of nonlinear equations in the setting of Banach spaces. First, 
we discuss the convergence analysis of our method by separating the equation into functions 
from which it is derived and the remaining part of the equation, according to the Generalized 
Theorem[1] for solving polynomial and transcendental equations. Solving an equation, either 
polynomial or transcendental, is solved only by categorizing the roots and not by some proce- 
dure of searching in individually intervals and this way was followed in the past by all existing 
methods. But the methods developed in the past are limited to isolated intervals without general 
acceptance. Finally, the method is strengthened by Newton’s method to speed up the finding of 
the roots. At the end of the analysis we give several examples with application of the method. 


Keywords: nonlinear equation, polynomial and transcendental equations, iterative Fixed Point 
Method, continuity conditions, Newton type method, initial values of roots. 


Part I. 
Introduction. ” Definitions & basically Theorems.” 


I.1. Inverse Function. The notion of an inverse is used for many types of mathematical construc- 
tions. For example, if we define f : T — S is a function restricted to a domain S and also a range T 
in which it is bijective and g : S + T is a function satisfying f(g(s)) = s for all s € S, then g is the 
unique function with this property, called inverse Function. It also follows that for all, so, i.e., inversion 
is symmetric. However, ”inverse functions” are also commonly defined for functions that are not bijective 
(most commonly for elementary functions in the complex plane, which are ), in which case, one of both of 
the properties may fail to hold (only for genuinely increase function). 


f (fF -'(y)) = y&y € f(T) 
Apply the forms { f-!(f(x)) = 2&a €T 


f (F-"(@)) = f@) 
Inverses are also defined for elements of groups, rings, and fields (the latter two of which can possess two 


different types of inverses known as additive and multiplicative inverses). Every definition of inverse is 
symmetric and returns the starting value when applied twice. 


I.2. Fundamental Algebra Theorem. 
If f(x) € F(x) and deg f(x) > 0, then there exists a € C such that f(a) = 0. That is, any non-constant 
polynomial with coefficients from the set of complex numbers have at least one root in C. 


I.3. Definition.I. If f(z) € F(x) where F' field, with deg f(x) = n, is analyzed in product of linear 
factors (splits) in F[x] if f(a) = c(a— a1) - (a — ag) +++ (a — an) where c, a1, d2,...dn € F. 


].4. Definition.II. Let E/F be a field extension and a € E. The a is called algebraic above F if f(x) € F(x) 
so that f(x) = 0 and f(@) = 0. If @ is not algebraic over F, then @ is called transcendental by F.An extension 
E of F such that the polynomial f(x), when we call it as the polynomial of E[z], has all its roots in E, i.e. 
analyzed in a product of a linear product to E/x]. When we talk about fields, E is a field extension of roots. 


I.5. How they are defined as Radicals & Periodic Radicals as Roots of f(x) = 0. If poly- 
nomial f(x) € C(x) with deg(x) < 5 has a field of roots E, it has as any root, a sum by finite count of 
radicals. If deg f(x) > 4 except for the cases analyzed in a product of linear factors, the roots take the form 
of sum of infinite radicals.Also as roots we have Sum of infinite radicals and for transcendental equations 
as we’ ll see above. 


1.6. Significant Definition III.” Categorization of roots in transcendental equations.” 

To find roots in transcedental equations we do not follow the procedure as in polynomials down to the 
fourth degree or Galois theory or other special method. We categorize each mononym function from 
which consists of the equation and after that we find the corresponding inverse function. 


The corresponding theorem is on page 5 and the formulas for the known functions on page 12. Of course 
and for polynomials up to 4th degree we use the same theorem. 


1.7. Cyclotomic Polynomials. Consider the cyclotomic polynomial «” — a = 0. We know that a 
polynomial is solvable if and only if its Galois group is solvable. There is a general method (we will explain 
in more detail below) for solving polynomial equations but other cases more complex with this solution of 
the equation. And this fact is because it involves the highest degree n. In the field of complex numbers, the 
solution of cyclotomic equation where n is an integer, are given by DeMoivre’s theorem: 


tp = Va: Ch = Va: e27**/” where k = 0,1,2,3...n—1. 


1.8. Iteration Method. In computational mathematics, an iterative method is a mathematical pro- 
cedure that uses an initial value to generate a sequence of improving approximate solutions for a class of 
problems, in which the n-th approximation is derived from the previous ones. A specific implementation 
of an iterative method, including the termination criteria, is an algorithm of the iterative method. An 
iterative method is called convergent if the corresponding sequence converges for given initial approxima- 
tions. A mathematically rigorous convergence analysis of an iterative method is usually performed; however, 
heuristic-based iterative methods are also common. Iterative methods are often the only choice for nonlinear 
equations. However, iterative methods are often useful even for linear problems involving many variables 
(sometimes of the order of millions), where direct methods would be prohibitively expensive (and in some 
cases impossible) even with the best available computing power. 


1.9. Attractive fixed points. 

If an equation can be put into the form f(x) = x, and a solution x is an attractive fixed point of the function 
f, then one may begin with a point x, in the basin of attraction of x, and let tn4, = f (xp) for n > 1, 
and the sequence {x,,},,., will converge to the solution x. Here x, is the nth approximation or iteration 
of x and x, 4, is the next or n+ 1 iteration of x. Alternately, superscripts in parentheses are often used in 
numerical methods, so as not to interfere with subscripts with other meanings. 


For example, x("+!) = f (2), If the function f is continuously differentiable, a sufficient condition 
for convergence is that the spectral radius of the derivative is strictly bounded by one in a neighborhood of 


the fixed point. If this condition holds at the fixed point, then a sufficiently small neighborhood (basin of 
attraction) must exist. 


I.10. Definition IV (Polynomials). Let F C C be a field and x element of C, algebraic above it 
F. We will say that element x is expressed with radicals if it belongs to a field E such that there is a 
sequence of fields 


F=FoCh,Choc...ckhCh4ic...ck ck 


for some integers s, where if we assume 2g an initial value of x we can take 
Fi41 = Fj (4/20) 


for some xo € Fj; and n is Natural number, 0 <i < s—1. With ~/zo we denote a root of polynomial 
x" — x9 = 0 € Fi [x]. The sequence is called a radical sequence of fields, and the E/F extension is called 
a radical extension. A polynomial f(x) € F'[] is called resolved by radicals above F, if any it is 
expressed in radically, that is, if any a radical extension containing the analysis of f(x) above F’. 


J.11. Galois theorem. Let F Cc C be a body, f(x) € Fla] with a field L. The polynomial f(x) is 
solvable above Q although only if the Gal (L/F) group is solvable. (At this point we have to give an 
explanation. We will not consider here the criteria for an extension L to be solved, so we’re not interested 
in the fact that there is a symmetric polynomial that consider Galois theorem. We will solve by repeating 
each set of roots in a general polynomial form. So an extension of our own is not necessary to be resolved 
radically by Galois theory. Therefore, our extension L’ by radicals may not be solved by Galois theory but 
by the method of repetition that we accept). 


Part II. 


I.1. Theorem 1.Generalized Theorem of existence and of finding Roots to a Equation.(GTRE) 
We assume that we have general, elementary functions (in general transcendental) that are analytic except 
for some isolated singularities and branch cuts, in which case these and their local inversions will have 
convergent Taylor series expansions on suitable disks. Under these conditions i consider it a transcendental 


f :C > C where f(x) = SO ajoi(x) + ao, ai,a0 € C with o;(x) functional factors where o; : C > C 
i=l 


for which there is likewise an inversion and if suppose G the field of the roots of equation f(x) = 0, that 
defined from relation by category of functional factors o;: 


G=F(G,= {ee ee G2 = {77 io | en) hes eee ate ss eee ck gece 
Sete st” ee cee 8) 


with 1 <<k<nand 0 < m, < o where G; subfields of total field of roots where produced the o;, and G 
it will be at the same time and the total global field, which belong the roots of the equation f(x) = 0. 


If we define the relation of function and divide it into 2 parts ie. 
1 n 
on(2) = up = — | - S- ajo;(x) — ao 


a 
k pT Ake 


then the subfields of roots that produces the o;(2), 1 < k < n (according to the generalized theorem by 


Lagrnage method and the (part#1(pag.1) that we mentioned [1]), and they will correspond to G,. Now 
if we accept all the functions o,(2) with 1 < k < n and exist inverse functions (invertible procedure) i.e 
are 1-1, the equation f(x) = 0 is fully resolved, and is represented by the field G. The counting of Roots 
of equation will be given from the relation |G] = |G, UG2U...UG,| where (Gi, G2,...Gn) the individual 
ones subfields) and will be a bounded integer or an infinite integer. 


Proof. 


n 
According to the equation S> aja; (x) + ao = 0,a;,a9 € C,i € NT apply the relations: 
i=1 


2 1=1,i42 
CE) coditae Sea vecs dha datentd mehagidvaer ener edeiotehanadanatnet 
1 n 
On(L) = UR = — [i- S- aioi(L) — ag 
ot i=1itk 
1 n-1 
On(L) = Un = — | - S- ajo; (x) — ao 
On i=1,ign 


From relations (I) and the property of inversion and with the iterative procedure applying from the previous 
ones relations we will takes 


1 ” = 
oth, = OL (2 (- Sa (xj) - os) =o7' (m1) EG, 


i>1 


1 n 
rh oy ee bee S- ajo; (a4) — ao =0, (uz) € Go 
a2 i=1,142 
OE Seer ery eto tr e tein sore oer eer eeter sete t est yar tary T or Metre ernst 
ef 4 “ 
dio =o" —{- S- a0; (£5) — ao =0;,' (ux) € G; 
Ok 1=1,i¢k 
1 n-1 
on he, eae (2 (- S- ajo; (xj) — os) = 0,1 (tn) € Gn 
w=1 


with yz the number of repeats ,1 << k <n, 0< mp < cw, 1 < q < my which means that the count of roots 


it can to be found in interval [0,00) and apply in generally for number 7, of the coefficient m;, (for any 
subfield G, number of roots) will be ym, € N+. The transformations of relation (II) they produce exactly 
the subfields G1, Ga,...,Gx,...Gn which generally they constitute the G. The genitive form finding of roots 


of f(x) = > ajoi(x) +. a9 = 0, {x,a;,a9 € C,i € N*} with the Periodic Radicals will be done with the gen- 
i=1 


eralized function ®; (uz) = 0, ' (ux) resulting from every relations in simple o;,(2) = uz => xy =O) (te) 
where o;(x) (functional factors of equation 1 < k <n) and each root will given from 


where the symbol # shows repetitive (periodic) procedure of values. For a number of repetitions 
we find each time the roots mx, “ the number of repeats, of category k from the form ¢/* x until some 
final value, with approach of error < 10*(—p), p is integer and p >> 1. The following stopping criteria 
is used for computer programs: 1.) |xj;41—2;| < ¢, 2.) |f(a;)| < ¢ where « = 10%(—p), p > 1, all 
Integers. Repeat procedures for all subfields of roots, are expressed computationally by the relations that 


result from (II) forms with programming language in the form (III). Analytically we will have them relations 


m 1 . 
oe a Pi (2 (— Sante a os) ) 


on Eo, = O2.u2 “a, s- ajo;i(#) — a0 , 
a2 i=1,i42 
n 
oo th, = Pnun ~ _ > ajioi(F#) — ao with Ym, € N*t,k € {1,2...n}, the number of re- 
i=1,idn 
peats 


Obviously if we have a limited number of functions which to be 1—1 and be inverses functions we will limit 
ourselves to the subfields they produce as many they can. The counting of Roots of equation will be given 
from the relation |G| = |G; UG2U...UG,| where G1, Go,...G, (the individual ones subfields) where for 
every one G_ corresponds an one inverse oj and the final count of roots of solutions of f(x) = 0 results 
from the sum: 


n ™; 


|G| = iy S- qi (Total number - solutions) 


i=1 j=l 


also qi the number of roots per m; category in the i position respectively (any category m; it has subcate- 
gories with count number from 1 to mj) and finally the 7m, is the number of subcategories (per i) category. 


II.2. Subfields, categories of subfields and number of roots per category [1] 


The f(z) = Sa; - o:(z) — ao = 0,{2,0:,0a0 € C,i € N*} it has roots Zo, where sg € Z,q € N 
i=1 


and s, is a multiple parameter with q specifying the number of categories within the corpus itself associated 
with either complex exponential or trigonometric functions. The Ly,, is the subfield concerns only this 


form, that is to say (o,). But the k-subfield itself can have q categories. In case we want to connect 
the subfields of the root in the complex plane and the general solution equation per category we will have 
the general relation: 


L = ao 1 = Qo 1 22 Qo 
abies og (= 08 (oy (2 = of (o,* (| 
q Qk a Ak Ak Ak 
n 
1 Be ao 1 
——of,|o,'| —-—- y 040; (Xo) (1*) 
& Oe SOS or Ae 


1<k<n,0<s<.co,qe€ Nt where of the complement of a,. 


Now, for the generalisation of cases, because this k takes values from to 1 +n, consequently the count 
of the basic subfields of roots also will be n, and consequently the field of total of the roots of the 
equation is L and will apply: 


n Uk 
fom) ERM BAe = || Tea th tig NUE RS ny, (2*) 
k=1 q=1 


The parameter u, takes different values and depends on whether a term function of the equation is trigono- 
metric or exponential or multiple exponential or polynomial. Example for trigonometric is 1 < ug < 2 and 
exponential with multiplicity 1 degree is ug = 1. We will look at these specifically in examples below. s, 
is a parameter with q specifying the number of categories within the body itself that are associated with 
either exponential or trigonometric functions or polynomial. ” But what is enormous interest about formula 
(1*) is that it gives us all the roots for each term of the equation o(z) = 0, once it has been analyzed and 
categorized using the inverse function technique o,, t (xo, Sq) after we do the analysis separately in each case. 


II.3. Number of Roots by use sets. In its general form, the principle of inclusion-exclusion states [2] 
that for finite sets A1,...,A,, one has the identity and just let A; be the subset of elements of S. Then we 
have 


JA, U A2U-+-UAnl= So |Al— So [Aim A5| 


1l<i<n 1<i<j<n 
+ So |Agn.Ayn Ag] —+ +++ (-1)"*1 Ar 9 42+ ++ An 
1l<i<j<k<n 
JA U Ag UU An| = SO(-1)*** SY [yer Ay 
k=1 JIC[n]:|J|=k 


Here we have the equality G; = A;Vi € {1,2...n} and G = S. Clearly the subfields can be infinite number 
elements. The G'; are subsets of the G' set of Roots of the equation f(x) = 0. 
Number intersections. The proof is done by induction: 


Cl >) |All} =dan, Cl © 14: 4g] ] =2 =n- (wt 1)/2,...,C (AN Asn...) =R=1 


l<i<n l<i<n 


The total count = 2” — 1. 


Example: The principle is more clearly seen in the case of three sets (n = 3), which for the sets A, B and 
C is given by 


S=|AUBUC|=|A]4+|B] +|C]-|AN B|-|ANC|—-|BnNC|+|ANBnC| 


a 


Figure 1 


This formula can be verified by counting how many times each region in the Venn diagram (fig.1) figure is 
included in the right-hand side of the formula. In this case, when removing the contributions of overcounted 
elements, the number of elements in the mutual intersection of the three sets has been subtracted too often, 
so must be added back in to get the correct total. This theory applies in our case both to the union 
of the principal fields of the basic functions and to the subfields of the partial functions. 


Part III. Method Infinity Periodicity of Approximate Roots (GRIM) 
III.1. Introduction. 


This method applies to polynomials with real or complex coefficients as well as transcedental polynomial 
formats. For such polynomials, we know that all complex roots appear as coupled pairs. The method 
attempts to locate the zeros of such polynomials by searching for zeros in a repetitive process after being 
split into two segments, and then in one section we calculate the inverse function that is the polynomial 
from mononym to until a trinomial. Therefore, we make the following reasoning, based on the theorem of 
continuity and if a function it has derivative and inverse. 


From special to general method: Suppose we have the function f = Yoo = 0, (1) we di- 


vide the equation with a sum of two functions equal to zero. So we will have a ‘apace and a polynomial. 
This namely f1(2) — fo(~) = 0. Because we are looking to find x, it will be the same root for both functions 
that as a sum are an equation we can equalize the one function fj(~) = u and fo(a) = u then with the 
relation xz, = f; (wu) = f, ’ (fo(am)) and same time apply x, = fy '(u) = fo’ (fi (@m)) for kz count 
of Roots for the category of function fz. Finally, with periodic iterative radicals we will have the final 
relation at the root of the initial equation (1) where 2, = fe Ci @e)) r= = 1,2,...N for r repetitions 


of maximum number N and certainly in this case we call relation f; (a -> ax" , fo(x) = —anx”. 


According to the Theorem and method G.R.L-N[1] the total of the roots will be the Union of 
the roots that will be found by the solution for the method (GRIM) for each monomial of polynomial 


So aja" = 0. Therefore, because i have n monomials will be n Roots. It is of course equivalent to the 


solution of the monomial with the maximum exponent, that is n. Of course this does not apply to tran- 
scendental equations in general, because there we have another categorization with complex parameters. 


III.2. Theorem 2. Generalized Iterative Method approximation of Roots (GRIM) 


We assume that we have general, elementary functions(in general transcendental) that are analytic ex- 
cept for some isolated singularities and branch cuts, in which case these and their local inversions will have 
convergent Taylor series expansions on suitable disks. Under these conditions for the determination of the 


roots any equation f(x) = > ajo;(x) — ap = 0, {x,a;,a9 € C,i € NT} (1) where f : C > C which consists 
i=1 


from distinct monomials functions o;(2),1 <i <A a; € C and for which there is likewise an inversion, 
after we accept oo(x) = 1, exist a repeating relation for every subfield of roots Gi, Go,...,Gr,...Gn and k 
the counting of the roots for any G,,r € {1,2,...,n}. Now if we now assume that that 2 initial value for 
x then and from GRE it will hold more closely with indices and for the any roots of f(x) will apply: 


1 ag 1 “ 
——o% | o>" re ee Sey S- 040; (Xo) jl<r<nl<k<o 


where of the complement of o;. 


where ys the number of repetitions and is based on the iterative procedure of generalized theorem GRE. 
This relation is also the generalized recurrent periodically algorithm of finding of the roots of any equation. 


Proof 


For proof we go from the special to the general case. Suppose we have the function f [2] -y Qjo;(x)-—a9 = 


0 (1) where f : C > C,a; € C where a; constants. I separate the initial equation wh a sum of func- 
tions and a monomial function and then therefore the two parts they will have sum equal to zero. This 


means as example: fi(x) + fo(x) — ao = 0 where: fi(x) = aioi(x) A fo(a = Saute > 241 = 


fe (2-- -— vf ))) (3) which results from the inverse of the first function with the simultaneous repet- 
a1 


itive Sere The result gives one of roots of subfield G,. Because we have a iterative process, it will 
come up with iterations for the given root type: 


In general we have for m, solutions corresponding to an equation of the of o, function i.e: 


= Qo 1 = Qo 1 = 0 
Reg ate oe ee pe ES ie a Pe -..n 
Or hy, ~ Wr & =f (ft Gi, aN NG 


where fr the complement of f,. 


Here zo is initial value from x and which results after from jp repetitions, take into account the same 
and the other fields Gi, Go,...,Gy,...Gn. If I call of(x) = > aj - oi(x) is complement of o,(x) where 
i=1,iA¢r 

Xo initial value for x then and from GRrE (Theorem 1) will be apply for the any roots of f(x): 


from (1) 
n 1 n 
0.01 (2) = — pe ajo;(2) — os)  o1(x) = — (2 2» ajo;(2) — os) 
1 i< 
Se rL=0; — S> aici (x) — ao (5) 
ay ~ 
1=2 
1 n 
and then apply for repeat with simple procedure ,,.7j41 =0,' | -— S- Q4 + 0; (2;) — Qo (6) 
ar m 
iArji=l1 


We can now easily deduce the generalized relationship 


ifk=1,0<s<o,qeN* 


n 
Apply : apox(x) + 0% (2) — a9 =O A 0; (2) = » a,o;(%) where ox the complement of ox. 
i=1itk 
where of the complement of o;. 


where ys the number of repetitions, m, emerges from the number produced by categorizing the roots for each 


G,,r € {1,2,...,n}, field and xo initial value for x that is based on the iterative procedure of generalized 
theorem GRrE. 


Part IV. Solved examples with the method (GRIM) [2,3,4,5] 
IV.1.1. Solving of general transcendental trinomial 


a2" a? +aq-2-b7" —ay =0 (ai, 7i;G,0,0 € Chi=1,2A a9 €C) 


2 
IV.1.2. We have the general equation with relation f(z) = S> a;-0;(z)—ao = 0, {ri, G, ai, ao € Cyt = 1, 2}. 
1=1 


Of course in this case the following will apply: o;(z) = 2"! -a*"" A a9(z) = 2" -b*"". But because we have 
identical monomials we can simply consider the one case it is enough to choose one of the 2 functions, in 
this case the 1 nd function i.e o1(z) = 2" - a7" ,r1,q, € C. According to the theory, we can to map the 
first term function to a variable, suppose u and will apply: 


o1(z) = 2" -a®" =uerlogz+2™ loga = logu+t 2knl,q, 40,kK EZ (1) 
From (1) = (2) -log 2 + z@ loga = logu + 2kaI(2), I: the imaginary unit. 


If I call z% = y (3) then (2) becomes 


71 


The point here is to relate variable y to a known solvable function as equation. This function is the 
W-function. If we now call the relationship 


W (ki, 8) 


ty-e%=s (5) Sty=W(ki,s)Sy= : 


(6), {s,¢€ C,t £0, ki € Z} 
If I now logarithmize relationship (5), I obtain the relationship 


ty oY = 5 « logy + ty = log = + 2knI (7) 


Combining (4,7) we arrive at the relationship 


By combining relations (3,6,8,9) we obtain the final solution which will be of the form 


1. 
l a i Qhon-l 
rs Si -W ky, cee un Pe n (10), {ko,k1 € Z},kg =1+q-1 
qi: loga —_ 


11 


In the same way we calculate for the second kernel the relation 


a1 
ro logb ag \ | 2kawt : 
gia Cod ea | 


q2 


Up to this point we have developed the concept of the iteration process. But the basic GRE theorem shows 


that to complete the solution we will have 2 formations for each case separately. i.e 01(z) = 2"!-a*"" Aoa(z) = 
272. be? 


According to the theory GRrE we developed before we consider 2 transformations. Here we have two 
functions 01(z),02(z) ie o1(z) = 2" + a*" A o9(z) = 27 -b*”. For each case we need to find the inverse 


10 


function separately per function. Therefore we expect to have 2 subfields of roots £1, Lz and therefore the 
total solution of the equation will be L = Ly U Lz. [1] 


IV.1.3. Finding the L, subfield. 


The first roots of sub-field LZ; results from inverse function and we give for kernel the relation 


eile, 
l a 10 Qson-I 
aie) = 09 fu =| (HE) awa #2 it) -e a , {52,51 € Z},52 =0+qm-1 (10) 


qm: loga n 


So, using this method (GRIP) we will have 1 relation to find the set of solutions of the first subfield. We 
talk about qi Categories of infinite (by parameter s;) complex roots. 


IV.1.4. Finding the Lz subfield. 


The second roots of sub-field Lz results from inverse function and we give for kernel the relation 


a1 
ro logb a2 \ | spn : 
cate) =u tty =| (a ites) W (oe | "e€ 9 , {85,81 € Z},85 =0+q@-—1 (11) 


So, using this method (GRIP) we will have 1 relation to find the set of solutions of the second subfield. We 
talk about q2 Categories of infinite (by parameter s/) complex roots. 


IV.1.5. Iterative Method approximation of Roots. 


In our case we have 2 functions, with one Sub-fields per function, which means that for the set of so- 
lutions of the field L of the equation we have L = L; U Ly. The iterative process that applies according to 
the theory we have developed will therefore be summarised in the formula 


i) Roots set for for the sub-field LZ, that generated by the function o1(z) 
1 Qa 1 a 
Lig cabs BO 2 he fd OO ck ae fy OU, os de 
O1 Sq 1 (= Qy 1 1 a4 Q4 1 1 a4 


2 
1 LH Qo 1 
eset zl pee Oi 1<k<2,0<s<00,qeEN’, 
at (a"(S-2-denes))))))) —e 


Apply: a101(z) + of(z) — a0 = OA of(z) = D7 agai(z) =ag-2”- be” 


mice 
] a 41 Qson-1 
oy *(z) = —— ‘W | 51, ras -e% 5 {82,81 € Z},82=0+q—1 
qj: loga oi 


where of the complement of 01 by the coefficient of az. 


ii) Roots set for for the sub-field Lz that generated by the function o2(z) 


1 
log b a2 q2 Qshn-l 
a= | (<5) w (a Me 8) -e2 , {82,81 € Z},82 =O0+q—1 


where o§ the complement of a2 by the coefficient of a. 

IV.1.6. We will look at a real example. We will look at areal example. Solve equation rae +4a? = 2a. 
According to the theory developed above we will have 

i) Roots set for for the sub-field LZ, that generated by the function oj (z) 


We have the constants: 


Pike oe COs AT at af Ow Ee fe (O08 Ss ce, 
o z O71 O71 O71 O71 O71 ees 
ne Qa, Qa, Qa, Qa, Qa, 


2 
Apply: a101(z) + of(z) — ao = OA of (z) = D7 agai(z) =ag-2”- Be 


1 
] a an Qsgn-l 
a, '(z)= 2 ‘WI 81, OE e ynt .e 4 {52,51 € Z},8. =0+q-1 
qi: loga a 


where of the complement of 01 by the coefficient of az. 


ii) Roots set for for the sub-field Lz that generated by the function o2(z) 


We have the constants: 


1 
1 1 
-—o§ (os (22- £-¥-aim) )))))) k =2, 0< 8 <00,q ENF, 
a2 a2 a2 i=] 


1 
Apply: az02(z) + 0$(z) — a9 = 0A 0$(z) = Di anoi(z) = a1 - 2" - be 


Zit 
_ r9 log b a2 a2 Qsh aI ; 
a= | (aes) w(x ag wt) ‘e @ , {82,81 € Z},82 =0+q-1 


where o§ the complement of a2 by the coefficient of a. 
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To get tangible results, we need to give a value for a. So we randomly choose a = 5. Therefore with the 
theory, we will therefore construct 2 programs depending on each of the assemblies we are considering. 


I) Program 1. Calculation of Roots set for the sub-field L; that generated by the function oj (z) 


Constants: a := 5;nl := 1;11 := —1;n2 := —-1;12:=1;b:=ajal :=1;k2 := 1; a0 := 2a; a2 := 1; 
Clear [xr,y];Off[FindRoot::Istol,Divide::infy,Infinity::indet ,FindRoot::jsing]; 

h2 [x_] = (—a2*xn2* b’ (x*(12)) + a0) /al 

h [x_] := al*x’n1*a‘ (x*(11)) + a2*x A n2* b* (x*(12)) — a0; 

For [kl = —20,k1 <= 20,k1 + +, 


flu_]:= 

N ( ProductLog [k1,n1/11 * log[a]*u*(n1/11)]* ((11/n1)/logfal)) (1/n1)* Exp [2*k2*I*7/n1]] ; 
xr = N[Nest[f[h2[#4]]&, 1/100, 5]]; 

FQ2 = N{xr, 15] 


Print ["2(’,k1,”) =”,s = N[y]/. FindRoot [h[y] == 0, {y, FQ2}, 
WorkingPrecision—>20],”,”,” Approach = ”, N{h{s], 2]]] 


Results 1: 


x(—10) = 4.56016142960097 + 39.0281378491641 I 
x(—9) = 4.43082289116071 + 35.1274487203240 I 


x(—8) = 4.28673006208215 + 31.2281347715740 I 
x(—7) = 4.12419318593592 + 27.3309211920972 I 
x(—6) = 3.93801587719292 + 23.4370540435161 I 


; 
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(—4) = 3.46061196593289 + 15.6705560300157 I 
(—3) = 3.14039657964064 + 11.8113468918753 I 
(—2) = 2.73364128469454 + 7.98932994522154 I 
x(—1) = 2.20125675516346 + 4.22307532327606 I 


; 


; 


) 
) 
) 
) 
(—5) = 3.72062664044545 + 19.5488002623928 I 
) 
) 
) 
) 


x(0) = 1.00000000000000 

x(1) = 2.20125675516346 — 4.22307532327606 I 
x(2) = 2.73364128469454 — 7.98932994522154 I 
x(3) = 3.14039657964064 — 11.8113468918753 I 
x(4) = 3.46061196593289 — 15.6705560300157 I 
x(5) = 3.72062664044545 — 19.5488002623928 I 
x(6) = 3.93801587719292 — 23.4370540435161 I 
x(7) = 4.12419318593592 — 27.3309211920972 I 
x(8) = 4.28673006208215 — 31.2281347715740 I 
x(9) = 4.43082289116071 — 35.1274487203240 I 


x(10) = 4.56016142960097 — 39.0281378491641 I 


20 complex roots (per 10 conjugates) in absolute ascending order with approximation 10*(—20). 


II. Program 2. Calculation of Roots set for for the sub-field Lz that generated by the function o2(z) 


Constants: a := 5;nl := 1;11 := —1;n2 := —1;12:=1;b:=a;al :=1;k2 := 0; a0 := 2a; a2 := 1; 
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Clear|xr,y];Off[FindRoot::lstol,Divide::infy, Infinity: :indet,FindRoot::jsing]; 

h2 [x_] := (—al*x“nl*a‘ (x*(11)) + a0) /a2; 

h[x] := al*x/n1*a/ (x*(11)) + a2*x A n2* b’ (x4(12)) — a0; 

For[kl = —20,k1 <= 20,k1 +4, 

flu] : = 

N | (ProductLog [k1, n2/12* log[b]*u’ (n2/12)]* ((12/n2)/log{b]))" (1/n2)* Exp [2*k2*I*2/n2]|; 
xr = N[Nest/f[h2[#]]&, 100, 7]] 


FQ2 = N{xr, 15]; 

Print ["2(’,k1,”) =”,s = N[y]/.FindRoot[h[y] == 0, {y, FQ2}, 
WorkingPrecision— 20],”,” ” Approach = ”, N[A[s], 2]]] 
Results 2 


x(9) = 0.00353456403606972 + 0.0280219318116643 I 
x(8) = 0.00431446145318829 + 0.0314301534678917 I 
x(7) = 0.00539824183640218 + 0.0357740085279072 I 
x(6) = 0.00697236486206852 + 0.0414959454659245 I 
x(5) = 0.00939555374766368 + 0.0493658249852404 I 
x(4) = 0.0134370696673487 + 0.0608465656289249 I 
x(3) = 0.0210242783146740 + 0.0790744219809211 I 
x(2) = 0.0383388136452226 + 0.112048875481559 I 
ax(1) = 1.00000005670143 + 3.53918967781277*10-° I 
x(0) = 1.0000000000000 


) 

) 

) 

) = 0.0134370696673487 — 0.0608465656289249 I 
—5) = 0.00939555374766368 — 0.0493658249852404 I 

) = 0.00697236486206852 — 0.04149594546592451 

) = 0.00539824183640218 — 0.0357740085279072 I 

) = 0.00431446145318829 — 0.0314301534678917 I 

) = 0.00353456403606972 — 0.0280219318116643 I 
x(—10) = 0.00295349037853237 — 0.0252774449784375 I 


Also 20 complex roots (per 10 conjugates) in absolute ascending order with approximation 10’(—20). 


As it is known we have only 2 infinite categories per root case. By combining the values between the 
2 approximations and the initial value, we achieve the best possible refinement of the roots. Usually we 
have 2 categories of initial values. Those that are in the interval [0,1) and those that belong to the interval 
(1, infinity). 


IV.2. Solving of polynomial equation. 


n 
For the determination of the roots of polynomial equation f(x) = S> ajoi(x) — ag = 0, {%,a;,a9 € C,i € 
i=l 
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Nt}, which consists from distinct monomials functions o;(x),1 <i<nAa; € C, exist a repeating relation 
for every subfield of roots G1, G2,...,G,,...G, and n the counting of the roots for any G,,r € {1,2,...,n}. 


If we dissociate the initial function f(z) = 3 aioi(z) — a9 = 0 if we dissociate the initial function 
n—-1 ae Fi n : ao 
f(x) = Gnon(x) + > a;0;(x) — a9 = 0 to 2 parts then we have the form aj ,) = —7- a iti 44+ =< > 
= i ier A= 
i Sere (-2 es aya’, + 2) , {ar € C,k = 0,1,2,.n—1,1 <r <n €N*} according to 
Oar if¢rixl ar 


the iterative procedure. The process starts with r =n, field {G,,} and continues with r=n—1,r=n-—2 
fields G,_1,Gp_9... ete until they arrive all the roots of equation and the field G. 


Example 
Calculate the roots of the equation a,2" — a,x — ap = 0, {an = 1,41, a9 € C} with the method GRrE. 


According to the iterative procedure, apply: 


1 1 1 1 1/n 

n i ao Qhri i ao 

(n&j41) = — S ati yy + > pBjq1 = errr. | — y ar, + — , 
i ae an ee An 


{a;, =1,01;09 € C,k = 0,1, 2, .n'— 1,n€ NT} 
with method GRrE a simple program will be ,if initial value is x0 and for 4 repeated times: 
A te | oe is 
f2 [x_] = al*x + a0; 
h2 [x_] := f2[x]; 
h{[x] := fl[x] — f2[x]; 
f [u_] := uA (1/n)* Exp [2 « k*a*I/n] ; 
xn = Nest/[f[h2[#]]&, x0, 4] 


The answer in general form is: 


1/n 
L/n 1/n 
L4= etki /n (x + aye2™*l/ (« + aye2™htin (a0 geen (ao + aiao)'/") ) ites 


with k=0,1,...n-—1 


The method can, of course, be combined with the Newton method because it is local, yielding results 
with a rather very good approach. Complete examples of such can be seen at the end of the work. 


IV.3. Equations of category with Transcendental or Irrational roots to polynomial form. 


We look at the form 2” — px* + q = 0,{(m,s) € R,(p,q) € C) (2) again in case c, but now with an 
extension of the set Z to R of G.R.P-N. The process is similar to b, we just have sub-cases that depend with 
the number of roots in relation to the maximum exponent m of the equation. We will see the whole process 
with examples in particular. The difference is to the fact that the number of roots in the transcendental is 
the integer part of the Real maximum exponent with a value +1, the final value of which depends on the 
sign of the last fixed term q, of the equation. So if it is a sign of q < 0 then the count of roots x of 
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transcendental equation (2) will be x = IntegerPart|/max{m, s}] +1, different if q > 0 then x = IntegerPart 
{max{m,s}]. This fact appears to be clear in all cases where we can meet. And because we are looking 
always the binomial | 7% = u, if a = max{m, s} that like if ¥ =u>a=e « u!/*,u€ C we distinguish 


two basic cases: 


I) [a] = IntegerPart|[a] = 2n+1,n€Z 


i) u>Othen & =0,+1,...+ ({a] — 1)/2 and 
ii) u<O then k =0,+1,...,+({a] — 1)/2,—([a] + 1)/2 


eR 


II) [a] = IntegerPant [a] = 2n,n € Z 


iii) u>0 then & =0,41,...+4 ([a])/2 and iv) u< 0 then & = 0,41,...,—([a])/2 


Examples 

1) [a] = IntegerPart[a] = 2n+1,n € Z 

i)u>0 

The binomial z*” = 2 it has complex roots. This equation is a transcendental one and accepts the fol- 


lowing solutions: Apply [7] = 3 therefore k = 0,+1 because k < ({a] — 1)/2 = 1. These solutions are given 
more analytically by the relations below: 


zo = 2/69 = 1.24686 
zy = 2'/™e2tt/™ — 0.51888 + 1.337744 
2-1 = 2% e—27t/* — 0.51888 — 1.33774: 


ii)u<0 


The binomial z7 = —2 it has only complex roots. This equation has solutions: Because apply [z] = 3 
therefore k = 0,+1,—2. We accept more analytically by the relations below: 


zo = (—2)!/7e° = 0.67368 + 1.04920i, z_1 = (—2)'/"e727/™ = 0.67368 + 1.04920: 
zy = (—2)/7 27/7 — — 1.93439 + 0.175951, 29 = (—2)'/7e7?- 274/™ — 1.23439 — 0.17595i 
We have 4 solutions 
II) [a] = IntegerPart[a] = 2n,n € Z 
iyu>0 
The binomial z™t! = 2 we ask complex roots. We accept the following solutions: Apply [7 + 1] = 4 


therefore k = 0,41,142 because k < ([a])/2 = 2. These solutions are given more analytically by the 
relations below: 


zg = 24D 6° = 1.18218 
2-1 = QV tV #27" +1) — 0,06345546 + 1.18047843i, zo, = 24/™e*47*/™ — —1.17537055 + 0.12672797i 


We have 5 solutions. 
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ii)u<0 


The binomial z™t! = —2 it has only complex roots. We get the following solutions: Apply [m + 1] = 4 
therefore k = 0,+1,—2. These solutions are given more analytically: 


zo = 21/(* +) 6° — 0.85807105 + 0.81318508i 

zy = QM (r+) e—2nt/(n+1) — 085807105 + 0.81318508i, 

zy = 21/7 6?™t/™ — _ 076595450 + 0.90048299i 

zg = 2/Te—47t/™ — _0.76595450 — 0.90048299i 
Therefore we have 4 solutions. 


At the same time, we also have the solution of the principal transcendental equation we shall need subse- 
quently for solving the more general forms of transcendental equations. 


IV.4. Equation category to full trinomial of the form a* + p-b*-—q=0,{a41Ab41,a,b,p,q € C}. 


We solve the two transcendental binomials a” = u and b” = uw that arise if we solve each term indi- 
vidually. The resulting fields for each binomial therefore arise are.. x = (logu+2-k-a-I)/loga and 
x = (logu+2-k-7-TI)/logb. 


We have 2 Programs in mathematica, for equation of form a? +p-b* —q=0,{a41Ab¥41,a,b,p,q € C} 
of G.R.P-N Method. Will be: 


1). 

Constants: a: = 7/5;b := 4/5; vu := 20;k :=0,+1,—-1+ 2,-2... infinity 
p:=Lq:= 1; 

hl x_]:=(a) 2+ p*(b)*x—4 

h2[x_]:=—p*(b)x + q 


k := 1;f [u_] := (Log[u] + (2*k**I)) /Log[a] 

xr = Nest/f[(h2[#])]&, 1, v]; 

s2 = N[Abs[xr], 5] 

FQ2 = N{xr, 10] 

Print|” Approach=” ,f1[d]/.FindRoot| [f1[d] = 0, {d, fq}, 
WoringPrecision— 10],”,2(”,k,”) =”, N[d]/ FindRoot [fl[d] = 0, 
{d, fq}}, WorkingPrecision— 10]]// Timing 


According to the program and these values for k i get: 


Results 3: 


k = 0,x = —0.444738772955 + 3.753482934416i 
k =1,x = 1.3939724221 + 17.6753885427i 

k = —1,x = 1.3939724221 — 17.6753885427i 

k = 2,x = 1.3685154337 + 38.41134620221 

k = —2,x = 1.8685154337 — 38.4113462022i 
Infinity Roots 


Significance has all the roots and which, as we know, is infinite. To complete the set of roots we should 
look at the other group of roots corresponding to the other field. 
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2). 


Constants: a := 7/5;b := 4/5; uv := 20;k :=0,+1,—-14 2, —-2...infinity, p:= 1;q := 1; 
hl [x_] = (a)*x-+p*( b)*x—q 

h2[x_] := —1/p*(a)x + q/p 

k := 1;f [u_] := (Log[u] + (2*k**I)) /Log[b 

xr = Nest [f[(h2[#])]&, 1, v]; 

s2 = N[Abs[xr], 5] 

FQ2 = N{xr, 12] 

Print|” Approach=” ,f1[d]/.FindRoot [fl[d] = 0, {d, fq}, 

WoringPrecision— 12],”,x(”,k,”) =”, N[ d]/.FindRoot [fi[d] = 0, 

{d, fq}, WorkingPrecision—12]]//Timing 


According to the program and these values for k i get: 


Results 4: 


k = 0,x = —1.89895932741 — 28.1072354396i 

k = 1,x = —1.89895932741 + 28.1072354396i 

k = 2,x = —0.528066150559 — 59.9081922009 
k = —2,x = —0.528066150559 — 59.9081922009 
Infinity Roots 


In roots calculation many times one of the two cases does not yield roots. Generally, however, the whole of 
the roots are the union of the roots of the two fields. 


IV.5. Solving the equation e® — 92? —5 = 0,x € C in the set of complexes. [7,11] 
According to the method of periodic Radicals for the solution we are looking for: 


1). To determine the functional factors 


In our case we have 2 functional factors o;(x) = e® A o2(x) = x?. Therefore the subfields that pro- 


duced are G1, G2 where G = G U G3 is the set of solutions of Roots of equation. 
2). To find each time the inverse function as to x, functional factor o,(x%),1<k <2 


For each functional factor resulting from the relation o,(x) = up = x9" = of (up) we take analyti- 
cally: 


2Qm1-r-i 


uy > ia; = (log(u)+2-w-m-i)/3-e°™ 3, 0< mm <2,we Z 


GO 0(t) = hes = Oa (ui) >e 


2 2Qm-n-t 


b) .09(@) Stig > 05,0; a5" (uz) > on" 25 = Un => OG a; = Une 2 ,O< me <1 


For both cases j denotes a repetitive process. 


n 
3. We determine approximate now in agreement with form bie atari = pu x (- x ajo; (Ft) — ao 
i=1i 


where the symbol # shows repetitive (periodic) procedure of values and 1 <= k <= 2 and y the number of 
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repeats. In our case we will have for each case separately: 
a) .°, Case: 


Qemy-rei 


i) Oy, (uy) = (log (uy) +2-w-m-d)-e ~F 0< mj S2,k=1wWEZ 
ii). #. ieee =, (2 (- de aioi(#) — «o)) =, (2 (—a202(#) — ao)) = ©, ,, (#? +5). 


1... Case: 


2-mg-n-i 


i) 1. Oo, (u2)=ftg-e 27 ,OSme<1,k=2. 


ii) pt er aaa = Din (2 (- > aj,oi(#) = «)) = O21, (2 (—a101(#) — ao)) = Do, (e# _ 5). 


i=1,142 


4. Numerical calculations for Complex roots: 

4.1. Categories of complex roots corresponding to the case a: 

According to (case a), there are 3 categories of roots and more specifically we take. 
The basics Programs in mathematica. 


d := 6; w := 0,1, 2; 

h2[x] := x2 +5; 

h [x_] := e*°3 — x42 —5; 

For [k = —d,k <=d,k +4, 

f [u_] := (2Imk + log [u])!/3 * Exp [2*I*1*w/3] ; 

xr = Nest/ffh2[#]]&, 1, 15); 

FQ2 = N{xr, 25]; 

Print|”x(”,k,”) =”, N[y]/ FindRoot [h[y] == 0, {y, FQ2}, WorkingPrecision — 20]]} 


i. If w = 1 for the first 13 roots (from infinite complex Roots) we get. 


Results 5: 

x(—6) = —0.0541597908772764481 + 3.4430687795713415686 I 
x(—5) = —0.0542272093684110848 + 3.2561465116970043623 I 
x(—4) = —0.0522799127977140002 + 3.0446951231985700472 I 
x(—3) = —0.044446198208043990301 + 2.7986415748311100303 I 
x(—2) = 0.073991657283343344493 + 2.2314465835450111411 I 
x(—1) = —0.043290439396162264200 + 1.8563014885856711080 I 


—0.51628371223967864689 + 1.0193275613354981293 I 
—1.4767721464906745539 + 1.0770106461916786832 I 
2) = —1.9196926480095895878 + 1.2656352955061933943 I 


) = 
1) = 
) = 
3) = —2.2255896653786903526 + 1.4135307520156373273 I 
) = 
) = 
) = 


x 
x 
4 
x(5 
x(6 


—2.4668779334538722208 + 1.5359648431918860639 I 
—2.6694198069847862912 + 1.6414453768855680311 I 
—2.8457158247398611933 + 1.7347795001319450124 I 


xX 


( 
( 
( 
( 
( 
( 
x(0 
( 
( 
( 
( 
( 
( 
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ii.If w = 0 for the first 13 roots (from infinite complex Roots) we get. 


Results 6: 


x(—6) = 2.9629015049765169380 — 1.6197988551044190385 I 
5) = 2.7970790204948353259 — 1.5157734611427453324 I 
4) = 2.6084406987474839689 — 1.3956240678529720229 I 
3) = 2.3870122571990767896 — 1.2512153777260514123 I 
2) = 2.1133652182429052922 — 1.0649772354315062827 I 


1) = 1.7402048569203715448 — 0.7832829915835906632 I 
= 1.2330961898317350649 


xX 


xX 


xX 


) 
) = 1.7402048569203715448 + 0.7832829915835906632 I 
2) = 2.1133652182429052922 + 1.0649772354315062827 I 
3) = 2.3870122571990767896 + 1.2512153777260514123 I 

i 

) 

)= 


x(1 


xX 


I 


x 
x(4 
x(5 
x(6 


2.6084406987474839689 + 1.3956240678529720229 I 
= 2.7970790204948353259 + 1.5157734611427453324 I 


oe 
ie 
= 
(— 
= 
(= 
x(0 
( 
( 
( 
( 
( 
( 2.9629015049765169380 + 1.6197988551044190385 I 


To this category belongs the only Real root of the equation. 


x(0) = 1.2330961898317350649 


iii. If w = 2 for the first 13 (from infinite complex Roots) we get. 


Results 7: 


6) = —2.8457158247398611933 — 1.7347795001319450124 I 
5) = —2.6694198069847862912 — 1.6414453768855680311 I 
4) = —2.4668779334538722208 — 1.5359648431918860639 I 

) 

) 

ion 


xX 


; 


3) = —2.2255896653786903526 — 1.4135307520156373273 I 
2) = —1.9196926480095895878 — 1.2656352955061933943 I 
1) = —1.4767721464906745539 — 1.0770106461916786832 I 
—0.51628371223967864689 — 1.0193275613354981293 I 
—0.043290439396162264200 — 1.8563014885856711080 I 
2) = 0.073991657283343344493 — 2.2314465835450111411 I 


a 
i 
VS 
3) = —0.044446198208043990301 — 2.7986415748311100303 I 
ye 
i 
pe 


xX 


x(4 
x(5 
x(6 


—0.0522799127977140002 — 3.0446951231985700472 I 
—0.0542272093684110848 — 3.2561465116970043623 I 


(= 
= 
(= 
oe 
(e 
Ge 
x(0 
( 
( 
( 
( 
( 
( —0.0541597908772764481 — 3.4430687795713415686 I 


4.2.Categories of complex roots corresponding to the case b: 
According to (case b), there is 2 Categories of roots and namely we have. 


The basics Programs in mathematica. 
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d:= 0;h2[x_] := e*3 —5;h[x_] := e*3 — x2 —5; 

For [k = —d,k <=d,k+4, 

f [u_] := u*(1/2)*e’ (2wzi/2); xr = Nest [f[h2[#]]&, xo, 6] ; FQ2 = N[xr, 25]; 

Print [” Approach =”, N[h[y],30]/ FindRoot [h[y] == 0, {y, #Q2}, WorkingPrecision > 30], 
”x(”,k,”) =”, Nly]/. FindRoot [h[y] == 0, {y, FQ2}, WorkingPrecision > 30]] 


If w = O&w = 1 we take 2 roots i.e we get 2 Complex Roots 


Results: 


x(0) = 0.073991657283343344493 — 2.2314465835450111411 I (with initial value xo = 1) 
&(0) = 0.073991657283343344493 + 2.2314465835450111411 I (with initial value xo = 1/2) 


The Approach in both cases is 10727 + 10774. 


IV.6. Solving the equation a®* + p-b” = q-c",{x,a,b,c} € C, in the set of complexes. In general, 
we have 1 methods of solutions and 2 cases. Both of these cases are generated using a logarithm. The 
functions that represent this method are cases that arise analytically per case: 


We define the function: 


Cc 


a) Case: for fi(a) =q—p-(8)* 


ae l wet 
i) Gay, (tn) = ETE w € Z. 


b) Case: for k = 2, fo(x) = ; _ 5 ; oe 


log wa+2-w-m-t 


i) ras O21, (u2) = log(2) 5 W Z. 


i) 05 F447 — Opp (2 (- Se eG aioi(H#) — a0) ) = ©, (fo(#)): 


Numerical calculations for Complex roots: 
6.1. Categories of complex roots corresponding to the case a: 


Constants: m, = 2 = W/43;m2 = 2 = — V3]; 


Cc Cc 


d:= 15;p:=7;q:=e; 

hi [x_] := (m1)*x + p*(m2)*x — q; 
h2[x_] := —p*(m2)x + q; 

For [k = 0,k <=d,k +4, 

f [u_] := (log[u] + (2*k**T)) /log[m1]; 
xr = Fold [f[h2[#])]&, —2, Range[10]]: 
FQ2=N [zar, 10]; 
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Print|” Approach=”, N{h1[y],30]/. FindRoot [hl[y] == 0, {y, FQ2}, WorkingPrecision— 30],”x(”,k,”) = 


”, N[z|/. FindRoot [h1[z] == 0, {z, FQ2}, WorkingPrecision— 30]]] 


For the first 15 complex roots (this case don’t have spouses complex roots) we get. 


Results 8: 

x(0) = 58.0487507536351199035328205089 + 7.61535990029372317354610795877 I 
x(1) = 1.32936292836325648633528741447 + 8.35263361944434710891163214802 I 
x(2) = 1.32936292837426009323673136176 + 16.7052672389407600393663302459 I 
x(3) = 1.32936292837426009323670176765 + 25.0579008584111400590493309734 I 
x(4) = 1.32936292837426009323670176765 + 33.4105344778815200787324412978 I 
x(5) = 1.32936292837426009323670176765 + 41.7631680973519000984155516223 I 
x(6) = 1.32936292837426009323670176765 + 50.1158017168222801180986619468 I 
x(7) = 1.32936292837426009323670176765 + 58.4684353362926601377817722712 I 
x(8) = 1.32936292837426009323670176765 + 66.8210689557630401574648825957 I 
x(9) = 1.32936292837426009323670176765 + 75.1737025752334201771479929202 I 
x(10) = 1.32936292837426009323670176765 + 83.5263361947038001968311032446 I 
x(11) = 1.32936292837426009323670176765 + 91.8789698141741802165142135691 I 
x(12) = 1.32936292837426009323670176765 + 100.231603433644560236197323894 I 
x(13) = 1.32936292837426009323670176765 + 108.584237053114940255880434218 I 


Infinity Roots 


6.2.Categories of complex roots corresponding to the case b: 
The basics Programs in mathematica. 


Constants: m; = $= V/43;m2 = b = —77/31; 
d:=15;p:=7;q:=e; 

hi [x_] := (m1)*x + p*(m2)*x — q; 

h2 [x_] := —1/p*(ml)*x + q/p; 

For [k = -d)k <=d,k +4, 

f [u_] := (log[u] + (2*k*a*T)) /log[m2]; 

xr = Nest [f[(h2[#])]&, 0, 10]; 

FQ2 = N{xr, 10]; 


Print|” Approach = ”, N[h1[y],30]/. FindRoot [hl[y] == 0, {y, FQ2}, WorkingPrecision— 30],”x(”,k,”) = 


”, N[z]/.FindRoot [h1[z] == 0, {z, FQ2}, WorkingPrecision — 30]]] 
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For the first 21 roots (from infinite complex Roots) we get. 


Results 9: 


10) = —17.6708670992525382500210144891 — 6.39244149043165467954365963797 I 
9) = —15.9052625797348216939204909555 — 5.74913079189852942523203026835 I 
8) = —14.1396581085817360423637483024 — 5.10582168159023406815737164189 I 
7) = —12.3740510102774718677920189344 — 4.46251796194974345001420000146 I 
6) = —10.6084256713969881208311412561 — 3.81922763811319342747954640732 I 
5) = —8.84271588555983753545329276412 — 3.17595005668124265151716787800 I 
4) = —7.07670147668304967743649663054 — 2.53256705974439390398062089406 I 
3) = —5.30984992287635303352948277895 — 1.88829785677614216529548120143 I 
2) = —3.54170969394503471525566086214 — 1.23955079654876832 102220047884 I 
1) = —1.77731279444547086200298016750 — 0.572568978830005695936055265763 I 


xX 


) = —0.0751212877365536830355891729551 + 0.152969515224423777278081181818 I 
x(1) = 1.32451209943781285438925570727 + 0.725825584824742783606616383329 I 
x(2) = 2.96729579227675451503974652116 + 0.814278965513704256780960100130 I 
x(3) = 4.89173812645201749253290815756 + 0.991485702541660584621184657797 I 
x(4) = 6.85122399379713247579236902591 + 1.22319558394859143399479007723 I 
x(5) = 8.81785390561844772357282670808 + 1.46634612658618689421058998067 I 
x(6) = 10.7864240033118422313395583168 + 1.71180405846075920210445681490 I 
x(7) = 12.7555184144055098977865402109 + 1.95770093244574709918109959409 I 
x(8) = 14.7247481449571241878636245713 + 2.20367427809457833024642873265 I 
x(9) = 16.6940113070916808934037897543 + 2.44965907041596002761446021192 I 


(ee 
(= 
( 
(= 
(= 
(= 
(= 
je 
e 
Ge 
x(0 
( 
( 
( 
( 
( 
( 
( 
( 
( 
( 


x(10) = 18.6632824183316702903725095336 + 2.69564502316182224459278506286 I 


Infinity Roots 


IV.7. Solving the equation ¢(x) = 2,a € C, in the set of complexes. [1,2,6,7,8] 


To solve the equation we use of the 2 functional equations of Riemann. We use 4 special algorithms per 
case. Riemann’s 2 functional equations are a function of ¢(), cos(), sin(), (27)”, Gamma( ) and depending 
on each function they produce the fields of the roots. 


C(x) — 2(2n)*—1 - sin (=) - Gamma(1 — 2) 


l-« 


C( 
C( 


Ta 
— 2(27)~* - cos (=) -Gamma(z) 
The functions that represent this method are cases resulting in detail per case: 


We define the functions (cases a & b) 


f(x) = Pee 2(27)*—* - sin (=) -Gamma(l— 2), ¢(4) =a 


C1 — 2) 2 
Case a: for k = 1, fi(x) = x/¢(1 — x) -1/2- (2-7)*(1 — x)/ Gamma(1 — x). 
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i)“. ®1, (ur) = (2 (a — arcsin (u,) + 2-7-w))/a,w € Z. 
si ee ae = Ory (2 (- det ajyoi(#) — ao) ) = 91, (fi(#))- 


i)“. ®2, (U2) = (2 (arcsin (ug) + 2-a-w))/7,w € Z. 
it) Ma = Gey (2 (-D4, aioi(#) — a0) ) = Ory (AH). 


I am using the function f(x) = 


Case b: for k = 1, fo(a) = ee ; er 
2 


We define the functions (cases a & b) 


bi) 
i) 7. Bay (uh) = 1/2 + (log (uf) +2-w- mT) /log(n), w € Z. 
i) mag = Oy (Fal): 


Case b: for k = 1, fo(x) = aie 5 
2 
bg) 
i)". Bay (uh) = 1/2 + (log (up) + 2+ w+ wT) /log(n), w € Z. 
ii) Mar) = 25, (fH). 


Numerical calculations for Complex roots: 
7.1. Category of complex roots corresponding to the case a(a1): 
The basics Programs in mathematica. 


Unprotect[G1]; 

fl[x_] := x/ Zeta[1 — x]*1/2*(2*7)*(1 — x)/ Gamma|1 — x]; 

f[z_] := x — 2* (2*7)* (x — 1)* Sin [1*x/2]* Gamma[1 — x]* Zeta[1 — x]; 

h[x_] := Zeta[x] — x; 
Off[FindRoot::Istol,FindRoot::cvmit,General::stop,FindRoot::bddir,FindRoot::nInum,ReplaceAll::reps, 
Infinity::indet ,FindRoot::srect,General::munfl]; 

Gl1[d_, q_,t_,s_,sl_,m_] := 

For [k = —d,k <=5,k+4, 

Fl[u_] := N[(2(m — ArcSin[u] + 27k))/z]; 

xr = Nest/F1(f1[#]]&, q, t]; 

FQ3 = N{xr, 80]; 

s3 = N[y]/.FindRoot [fly] == 0, {y, FQ3}, WorkingPrecision > m|] 

s4 = N[yl]/.FindRoot [f[y1] == 0, {y,s3}, WorkingPrecision > mJ; 

If [Abs[s3 — s4]] < s && Abs|s4] < s1, 

Print [” Approach=”, [h[y1], m]/. FindRoot [h[y1] == 0, {y1, s4}, WorkingPrecision > ml], 
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%a (",k,”) =”, N[y1]/.FindRoot|h[y1] == 0, {y, s4}, WorkingPrecision — ml] , Loopback] 
Protect[G1]; 
G1[20, 1/2, 10, 0.00001, 75, 70] 


Results 10: 


x(—9) = —33.999999999683437233174912281325947 
x(—8) = —29.999999800525511561320078724806606 
(—7) = —25.999927036563283322479328779219027 
(—6) = —21.985531995733483012809206667804255 
(—5) = —17.729166140171729124345484517377206 — 1.1113876772379532029008867533810019 I 
x(—4) = —14.613492351835987588071919186650501 + 3.1007918168141534354604138798653371 I 
(—3) = —10.821148006023595165258300978174459 — 5.3193414062340658907809379747379031 I 
(—2) = —5.2793950307843170303114858336796244 + 8.8027675740832589731484730202165154 I 
(—1) = —0.29590500557521395564723783108304803 
x(0) = —0.29590500557521395564723783108304803 
( 


x(1) = 1.8337726516802713962456485894415236 


Infinity Roots 


The approach from top to bottom is of the order of 10*(—22) to 10(—36). 


7.2. Category of complex roots corresponding to the case a2: 
The relevant basics Programs in mathematica. 


Unprotect[G1]; 
fl[x_] := x/Zeta[1 — x]*1/2*(2*7)*(1 — x) /Gamma|1 — x]; 

f[x_] := x — 2*(2*7)*(x — 1)* Sin [7*x/2]" Gamma[1 — x]* Zeta[1 — x]; 
h [x_] :=Zeta|x] — 

G1[d_,q_,t_,s_,sl_,m_] := 

For[k = —d,k <= 2,k4 


% 


F2[u_] := N[(2(ArcSin[u] + 27k))/7]; 
xr = Nest/F2[f1[#]]&, q, t]; 
FQ3 =N[xr,80]; 


s3=N[y]/.FindRoot|fly] == 0, {y, FQ3}, WorkingPrecision > m]; 

s4 = N[yl]/.FindRoot [f[y1] == 0, {y1l,s3}, WorkingPrecision—m]; 

If [Abs[s4] < s1,Print [’x(”, k,”)=”,N[y1]/.FindRoot [h[y1] == 0, {y1, s4}, 
WorkingPrecision—m]],Loopback]] 

Protect[G1]; 

G1[20, 1/2, 10, 0.00001, 75, 70] 
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Results 11: 


x 


8) = —32.000000008467671636378705967757058 
7) = —28.000004102805105685691855067848468 


x 


x(—6) = —24.001107180111899006547352622573465 

x(—5) = —20.131186308008286295709692618189034 

x(—4) = —14.613492351835987588071919186650501 + 3.1007918168141534354604138798653371 I 
x(—3) = —10.821148006023595165258300978174459 + 5.3193414062340658907809379747379031 I 
x(—2) = 1.8337726516802713962456485894415236 

x(—1) = —5.2793950307843170303114858336796244 + 8.8027675740832589731484730202165154 I 


( 
( 
( 
( 
( 
( 
( 
(—1) 
w(1) = —0.29590500557521395564723783108304803 
w(2) = —0.29590500557521395564723783108304803 


The approach from top to bottom is of the order of 10*(—22) to 10“(—36).In short, we have 4 complex 
in this area(categories)and the rest are real. 


7.3. Category of complex roots corresponding to the case bl: 
The relevant basics Programs in mathematica 


Unprotect [G1]; 

Clear All; 

f2[x_] := x/ Gamma|(1 — x)/2]*Gamma|x/2]/(Zeta[1 — x]); 

f[x_] := Gamma|x/2] — (Zeta[1 — x]/x)*a*(x — 1/2)* Gamma|(1 — x)/2]; 

h [x_] :=Zeta [x] — x; 

Off[FindRoot:: stol,FindRoot::cvmit,General::stop,FindRoot::bddir,FindRoot::nlnum, 
ReplaceAll::reps, Infinity: :indet,FindRoot::srect,Less::nord,Divide::infy,Set::setraw]; 
G1[d_,q_,t_,s_,sl_,m_] := 

For [k = —d,k <=d,k+4, 

F3 [u_] := 1/2 + (log[u] + 2*k**T) / log[z]; 

xr = Nest/F3[f2[#]]&, q, t]; 

FQ2 = N{xr, 10]; s3 = 0;s4 = 0; 

s3 = N[y4]/.FindRoot [fly4] == 0, {y4, FQ2}, WorkingPrecision— 30]; 

s4 = N[yl]/.FindRoot [f[y1] == 0, {yl, s3}, WorkingPrecision— 30]; 

If [Abs[s4] < sl, Print|” Approach =”, N{ hly],m]/ FindRoot [h[y] == 0, {y,s4}, 
WorkingPrecision— m],”x(”,k,”) =”, N[y]/.FindRoot [h[y] == 0, {y, s4}, 
WorkingPrecision— mJ], m]]; Loopback |]// Timing 

Protect[G1]; 


G1(20, 1/2, 10, 0.00001, 110, 30] 
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Results 12: 


x(—14) = —1.26139081555873261406489769808 — 87.7583685099246485627976563544 I 
x(—13) = —1.12616495962146742040656000337 — 80.4582743596358549469056379024 I 
x(—12) = —1.29290215735362288709093453220 — 75.6417409652417064465387628638 I 
x(—11) = —1.26227506725207853187875606381 — 70.3556950691515434161289031984 I 
x(—10) = —1.29443568561644493920364442108 — 65.2301042074151417066393823324 I 
x(—9) = —1.41821905056698762133825325680 — 59.6454164781995695989688716976 I 


=9) = 
—8) = 
-7)= 
—6) = 
—5) = —1.51767569589846950803032259371 — 38.1936599169588114867812150952 I 
Ay= 
f= 
ys 
ae 


8 


—1.46547485014965062526389943588 — 48.0489077335066889394342997498 I 


‘al 


8 


8 


4 
3 
2 
1 


8 


—1.6934543711272871415911807768 — 26.5282868639131487782281764409 I 


8 8 8 


LOS LL LON ON ESS ON FS EE ES FE OS IO IS ES EN ENO REIN Le a 


) = 1.83377265168027139624564858944 

x(1) = —1.74955435401611051766624316629 + 30.8158449441671134562051307077 I 
x(2) = —2.03689852362822026102533100727 + 21.993069828431773685 1003882855 I 
x(3) = —1.6934543711272871415911807768 + 26.5282868639131487782281764409 I 
x(4) = —1.74955435401611051766624316632 + 30.8158449441671134562051307077 I 
x(5) = —1.51767569589846950803032259371 + 38.1936599169588114867812150952 I 
x(6) = —1.55341651249244969842510775727 + 41.5258170026898503231006830279 I 
x(7) = —1.46547485014965062526389943588 + 48.0489077335066889394342997498 I 
x(8) = —1.28703179907482141334672081309 + 53.9110489785536948914586838211 I 
x(9) = —1.41821905056698762133825325680 + 59.6454164781995695989688716976 I 
x(10) = —1.29443568561644493920364442108 + 65.2301042074151417066393823324 I 
x(11) = —1.26227506725207853187875606381 + 70.3556950691515434161289031984 I 
x(12) = —1.29290215735362288709093453220 + 75.6417409652417064465387628638 I 
x(13) = —1.12616495962146742040656000337 + 80.4582743596358549469056379024 I 
x(14) = —1.26139081555873261406489769808 + 87.7583685099246485627976563544 I 


Infinity Roots 
7.4. Category of complex roots corresponding to the case b2: 


The relevant basics Programs in mathematica 


Unprotect[G1]; ClearAll, 

f2 [x_] := x/ Gamma|(1 — x)/2]* Gamma|x/2]/(Zeta[1 — x]); 

f[x_] := Gamma|x/2] — (Zeta[1 — x]/x)*a*(x — 1/2)* Gamma|(1 — x)/2]; 
h [x_] := Zeta[x] — x; 


—1.28703179907482141334672081309 — 53.9110489785536948914586838211 I 


—1.55341651249244969842510775727 — 41.5258170026898503231006830279 I 


—1.74955435401611051766624316632 — 30.8158449441671134562051307077 I 


—2.03689852362822026102533100727 — 21.9930698284317736851003882855 I 
—1.74955435401611051766624316629 — 30.8158449441671134562051307077 I 


Off [FindRoot::lstol,FindRoot::cvmit,General::stop,FindRoot::bddir,FindRoot::nlnum, 


ReplaceAll::reps, Infinity: :indet,FindRoot::srect,Less::nord,Divide::infy,Set::setraw]; 
G1 [d_, q_,t_,s_,sl_,m_] := 
For [k = —d,k <=d,k+4, 
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F4[u_] : 
xr = Nest 


—1/2 + (log[u] + 2*k**I) /log|[z]; 

A[f2[#]]&, q, t]; FQ2 = N[xr, 10]; s3 = 0;s4 = 0; 
s3 = N[y4]/.FindRoot [f[y4] == 0, {y4, FQ2}, WorkingPrecision — 30]; 
s4 = N[{yl]/.FindRoot [f[y1] == 0, {yl, s3}, WorkingPrecision— 30]; 
If [Abs[s4] < sl, Print/” Approach =”, N{ h[y], m]/. FindRoot [h[y] == 0, {y, s4}, 
WorkingPrecision— m],”x(”,k,”) =”,zl = Nly]/. FindRoot [h[y] == 0, {y, s4}, 
]; Loopback]]//Timing 


[F 
4] 
1 


WorkingPrecision— m 
Protect[G1]; 
G1[20, 1/2, 10, 0.00001, 110, 30] 


Results 13: 
xu(—14) = —1.23852071108706795556935851810 — 85.4090090247636352615544742032 I 
x(—13) = —1.29918130011957305015490858296 — 77.9865709309247937726779884803 I 
x(—12) = —1.15953302973918715196482277792 — 72.9767004265557497382345340656 I 
x(—11) = —1.33095139800651049656449887081 — 67.7834609908689492473936882734 I 
a(—10) = —1.1924742451865064471089288392 — 62.3193286218473504568063554479 I 
x(—9) = —1.3366144881440349550167977808 — 56.9002016278718488305446496687 I 
x(—8) = —1.40616347663699362699926390436 — 50.9143847707163058133311938748 I 
xu(—7) = —1.32728372617780452008783857748 — 44.7148157580246341928779393786 I 
x(—6) = —1.55341651249244969842510775728 — 41.5258170026898503231006830279 I 
x(—5) = —1.51383649144589060970515941062 — 34.4539606289612051139663453008 I 
u(—4) = —1.74955435401611051766624316641 — 30.8158449441671134562051307077 I 
x(—3) = —2.03689852362822026 102533100727 — 21.993069828431773685 1003882855 I 
x(—2) = —2.38593568342424543039818327574 — 16.2709870621967901723315117424 I 
x(—1) = —1.5138364914458906097051594106 — 34.4539606289612051139663453008 I 
x(0) = —1.55341651249244969842510775736 + 41.525817002689850323 1006830279 I 
x(1) = —1.5138364914458906097051594106 + 34.4539606289612051139663453008 I 
x(2) = —2.38593568342424543039818327574 + 16.2709870621967901723315117424 I 
x(3) = —2.03689852362822026 102533100727 + 21.993069828431773685 1003882855 I 
x(4) = —1.74955435401611051766624316641 + 30.8158449441671134562051307077 I 
x(5) = —1.51383649144589060970515941062 + 34.4539606289612051139663453008 I 
x(6) = —1.55341651249244969842510775728 + 41.525817002689850323 1006830279 I 
x(7) = —1.32728372617780452008783857748 + 44.7148157580246341928779393786 I 
x(8) = —1.40616347663699362699926390436 + 50.9143847707163058133311938748 I 
x(9) = —1.3366144881440349550167977808 + 56.9002016278718488305446496687 I 
x(10) = —1.1924742451865064471089288392 + 62.3193286218473504568063554479 I 
x(11) = —1.33095139800651049656449887081 + 67.7834609908689492473936882734 I 
x(12) = —1.15953302973918715196482277792 + 72.9767004265557497382345340656 I 
x(13) = —1.29918130011957305015490858296 + 77.9865709309247937726779884803 I 
x(14) = —1.23852071108706795556935851810 + 85.40900902476363526 15544742032 I 


Infinity Roots 


The approach from top to bottom is of the order of 10“(—25) to 10“(—30). 
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Epilogue 


As we can see, each solution of an equation must be organized by categorizing its roots, so that the full 
range of roots is covered and not individually at some intervals. It is not always easy to choose the initial 
value, but we know 2 regions the interval (—1,1) and the maximum value of the full function which is into 
region [1,r >> 1). The choice and arrangement of examples is due to colleague Etairi Biragova thanks. 
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